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Abstract 

The color superconductivity of a dense quark matter is reviewed with em- 
phasis on the long range nature of the pairing force and the multiplicity of 
the order parameter. The former gives rise to a non BCS behavior of the 
superconducting energy scale and the latter modifies the critical value of the 
Ginzburg-Landau parameter that separates the superconductivity of type I 
and that of type II. 
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I. INTRODUCTION 



In this lecture, I would like to review our works for the past years on the color supercon- 
ductivity in a dense quark matter. While it is an old subject dated back to 1970's [1] [2] 
[3], the revived interests on the color superconductivity since 1998 [4] [5] [6] are largely 
promoted by the experimental efforts of exploring the phase diagram of the strong interac- 
tion through RHIC or inside neutron stars. While lattice simulations have been successful at 
nonzero temperature, it becomes prohibitively difficult when the chemical potential becomes 
nonzero. Analytical technique remains the only means especially at large chemical poten- 
tial, where an ideal gas of quarks and gluons presents the leading order approximation. The 
important progresses that have been made recently include the discovery of the color-flavor 
locked structure in the super phase [6] and the first-principle determination of the energy 
gap and the transition temperature [7] [8] [9] [10] [11] [12] [13]. The entire subject has 
been reviewed by several authors [14]. 

Among many novel characters of the color superconductivity, I would like to take this 
opportunity to highlight two of them which differ remarkably from the non-relativistic su- 
perconductivity in a metal. The first is the long range nature of the pairing force, which 
is responsible to the non BCS exponent of the weak coupling formula of the gap energy or 
the transition temperature and the suppression of the pre-exponential factor. The second 
is the multiplicity of the order parameters, which leads to a rich variety of inhomogeneous 
condensates and a different value of the critical Ginzburg-Landau parameter that separates 
the superconductivity of type I from that of type II. Many technical details have been sup- 
pressed and this lecture will serve a tour guide to our papers on the subject and the related 
works by others. 
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II. THE LONG RANGE PAIRING FORCE. 



Let us consider a quark matter of ultra high baryon density such that the corresponding 
chemical potential is well above Aqcd to warrant a perturbative treatment because of the 
asymptotic freedom. While this density may not be attainable in a realistic quark matter, 
say the core of a neutron star, the approximation there is systematically under control. The 
Lagrangian density of QCD at a nonzero chemical potential // in the chiral limit reads: 

1 - d 

C — — trF^j^F^j^ — igAfj)ip + iJiil)lA''P + renormalization counter terms, (1) 

where A^, = A^T^ and F,,^ = ^ - ^ - i9[\, M with the SU{N^) generator in its 
fundamental representation and ^ is a Dirac spinor with both color and flavor indices. The 
renormalized coupling constant is defined at the chemical potential via 

9 = X (2) 

(117V,-2iV;)lnf 

with II » A and A = Aqcd for N^^ Nf ^ 3. 

Perturbatively, the di-quark interaction is dominated by the process of one- gluon ex- 
change, as is shown in Fig. 1. The amplitude is simply that of the one-photon exchange 
in QED multiplying the group theoretic factors jj'^i^i jj'=2'=2^ which can be decomposed into a 

color anti-symmetric channel (anti-triplet for SU{3)) and a color symmetric one (sextet for 
SU{3)), i.e. 

As both the electric and magnetic parts of the one-photon exchange are repulsive between 
two electrons fiying in opposite directions, the one-gluon exchange interaction between two 
quarks flying in opposite directions is attractive in the color anti-symmetric channel because 
of the negative sign of the first term on r. h. s. of (3). 
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FIG. 1. The one-gluon exchange vertex. 

At high baryon density, the Fermi sea of quarks tends to screen the one-gluon exchange 
interaction through the HDL(hard dense loop) resummed gluon propagator [15], which 
takes the form 



A;^ + a;^ + aM{k,uj) 



(4) 



(5) 



and D4j{'k,iu) = in Coulomb gauge with iu the Matsubara energy of the gluon. In the 
region of energy and momentum for pairing, i.e. k « fj,, cu « fj, and to « k, the magnetic 
self-energy is given by 



/7 \ ^ 2 \^\ 



(6) 



and the electric self-energy by 



aE{k,u;) c^m% (7) 

with m|) = the Debye mass. The absence of the screening of the long range mag- 

netic gluon propagator in the static limit, a; = 0, leads to the forward singularity of the 
di-quark scattering, which dominates the pairing force. The long range magnetic gluons 
introduces also the non Fermi liquid behavior, that renders the quark self energy function 
logarithmically enhanced toward the Fermi surface, i.e. 



-1 

-^,^_i.ln (8) 

with iv the Matsubara energy of the quark and Kc {ksT « Kc « n) cutoff momentum. 

The energy scale of the color superconductivity is measured by the energy gap of the 
quark spectrum at T = or the superconducting transition temperature Tg. The determi- 
nation of the latter quantity will be reviewed in this section. 

The proper vertex function r(iz/', p'|w, p) for di-quark scattering with incoming energy- 
momenta (p, and outgoing energy-momenta {p', satisfies a Dyson-Schwinger equa- 
tion, which, when projected into the color anti-symmetric channel with an angular momen- 
tum J, takes the form 

Vj{iu',p'\ip,p) = Tj{iu',p'\iu,p) + /c^T V / dqKj{iu' ,p'\w" ,q)Tj{w" ,q\w,p), (9) 
where rj(w',p'|w,p) stands for the two particle irreducible part of rj{ii'',p'\ii',p), the kernel 

Kj{iv ,p\iiy , g) = 2n^{2J+l) ^^'^^'P)^^~'^^'P> vlO) 
with S^iv, p) the full quark propagator. The spinor indices and their contraction structure 
have been suppressed in (9) and (10). To the lowest order in coupling constant, T is give by 
the one-gluon exchange diagram in Fig.l and the quark propagators in the kernel (10) are 
replaced by bare ones. For v « fj,, v', and p ~ p' ~ we have 

rA.^M - + ^) (1" ^ + 3c.) . -^(1 + -1) (i„ ^ + 3c.) 

(11) 

where z> = 2^(?) ^^^^ cj = 1 for J = and cj = exp ( - 2 En ^) for J 0. The last 
step of (11) follows from Son's approximation [7] with z/> = max(i/, u'). 

The pairing temperature within each angular momentum channel, Ty\ corresponds to 
the highest temperature at which the Fredholm determinant of eq.(9), Dj — det(l — Kj) 
vanishes and the transition temperature is Tc — max {T^^) for all J. Here we have assumed 
the pairing between two quarks of the same herlicity, which contains the s-wave channel 
{J — 0) and will take the maximum advantage of the pairing force. 



To highlight the impact of the long range pairing force in the weak coupling formula for 
Tci we employ the standard expansion of the Predholm determinant, 



^ Jo 2 Jo Jo 



K{ii',p\iu,p) K{ii',p\ih'' , q) 

+ 

K{ii>' ,q\ii>,p) K{iv' ,q\iv\q) 
(12) 



For pairing mediated by phonons, which is of a short range, the m-th term of the expansion 
is of the order of g^'^ In with ujd the Debye frequency. At the transition, In ~ 1 
and the remaining terms are of the order of g'^"^~'^. This gives rise to the scaling formula 
ksTc — cbcs exp(— ■^^^) with kbcs fixed by the term of m = 1 and the leading order of 
cbcs fixed by the term of m = 2. In case of QCD, the logarithm inside the kernel (10) 
makes the m-th term of the expansion go like g'^'^ In^"^ and the Frcdholm determinant is 
dominated by a function of gln^^. It follows then that the transition temperature scales 
hke ksTc — cqcd exp — ^^^^ [7] [16] [17] leaving both kqcd and cqcd non-trivial to 
determine. 

A perturbative method was developed in [11] [12], which yields the formula of the 
pairing temperature 



ksT^ = ccV'c,4e-f [1 + O(^ln^)], (13) 



where the exponent factor k — ^ jv^i "^^ ^^^^ obtained in [7], the pre-exponential 
factor c = 1024A^7r^7Vp was found in [8] and [9], the factor d = 2e^ with 7 the Euler 
constant was found in [9] and [12]. The factor c" = cxp — ]^(^^ + 4)(A''c — 1) stems from 
the non Fermi liquid behavior of the quark self energy (8) and was calculated in [12] and 
reproduced in [13] via the gap equation (The existence of this correction was suggested in 
[7]). The gauge invariance of the formula (13) was discussed in [18]. The angular momentum 
dependent factor cj was calculated in [12] and the consequence of the universal exponent 
for different angular momentum channels in (13) was discussed in [19] in the context of the 
crystalline superconductivity. Some higher order terms have also been identified [20]. It 
follows from (13) that the superconducting transition temperature is Tc = T^^K 
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A similar perturbative method for the gap energy has been developed in [21]. 



III. THE MULTIPLICITY OF THE ORDER PARAMETER. 

The Ginzburg-Landau theory is a powerful tool to explore the superconducting state right 
below the transition temperature. Unlike the ordinary superconductors, the order parameter 
of the color superconductivity which describes the di-quark condensate, '^J^f^ is of multi- 
components, carrying the color-flavor indices and the chirality. For three color and three 
flavors, the symmetry group of theory in the chiral hmit is SU{?>)c x SU{3)fj^ x SU{3)fj^ x 
U{1)b of which the electromagnetic gauge group, [/(l)cm is a subgroup. Restricting within 
the even parity sector and neglecting the small projection in the color-sextet representation, 
the order parameter takes the form 

= ^"'"'"^hhf^- (14) 

and the Ginzburg-Landau free energy functional which is consistent with the symmetry 
group reads [22] [23] 

J d^f[^FljFlj + xA)^ + 4tr(^$)^(^$) + 4atr$^$ + 6itr($^$)^ + fe2(tr$^$)^] 

(15) 

where the covariant derivative = V$ — igA^ + i-^eA^T^ with A'' the color vector 
potential and A the ordinary electromagnetic vector potential. The 3x3 matrix $ = 
$0 + ^iT'' with $0 the singlet and the octet under a simultaneous color-flavor rotation. 
The generator here, T^, should be understood to be that of the anti-triplet representation 
of SU{3)c. At weak coupling, the one-gluon exchange lends us the following expressions of 
the Ginzburg-Landau coefficients [22] [23] 



a — 



k%Tc{T-Tc), 



48n^ 
7<(3) 

2 



= (16) 
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which follows also from the mean field approximation of the NJL effective action at moderate 
coupling if ksTc « [x. 

Let us review first the structure of a homogeneous condensate [22], i.e. A = A = V$ — 0, 
for which three regions of the parameter space need to be considered respectively. 1) 61 > 
and hi + 862 > The minimum free energy corresponds to the color-flavor locked condensate 
[6] with $ = (poU , where 

and is a unitary matrix (The nontrivial winding of U onto the gauge group, S'C/(3)c x 
?7(l)cm. gives rise to vortex filaments [23]). 2) 61 < but 61 + 62 > 0: In this case the color- 
fiavor locked condensate (17) becomes a saddle point of the free energy, which is nevertheless 
bounded from below. The minimum corresponds to an isoscalar condensate [22], given by 
$ = diag(^^— ^^^^^ e'", 0, 0^. 3) For hi and 62 outside the region specified by 1) and 2), the 
free energy is no longer bounded from below. Higher powers of the order parameter have to 
be restored and the superconducting transition becomes the first order one. In the following 
we shall focus our attention to the case 1) of the parameters. 

The color-fiavor locked condensate breaks the original symmetry group down to a SU (3) 
of simultaneous color-fiavor rotation, among which there is an unbroken U{1) gauge sym- 

— * 

metry and the corresponding gauge potential V is obtained through the rotation [24] 

V = -^^sin e + Acos 9 

V =A^cose + Asm9, (18) 

9 — * — * 

with tan^ = In comparison with the electroweak theory, the field V and V are the 

analogs of the photon and the Z-boson, and the rest of the A^s are like the VF-boson's. The 
mass of the V and that of A' (Z = 1, 7) are 

2 A 2 i2 2n ^ 2 2 2 a ^ 

ruz = 0osec = ^, rriy^ = m^cos = 

with 5 and 5' the corresponding penetration depths. Other low-lying excitations of even 
parity consist of a Goldstone boson associated to the broken baryon number, [/(l)^, the 
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Higgs bosons associated to V and to A' (/ = 1, 7) with masses 

2 2 

= (^1 + 3^2)00 = ^> "^'1 = bi(f)l = 

and the corresponding coherence lengths. The excitations of odd parities includes the 

Goldstone bosons associated to the chiral symmetry breaking triggered by the color-flavor 
locking [6] . 

An interesting issue to address is the type of the superconductivity in response to an 
external magnetic field, which is partially screened because of its projection onto the broken 
V field through (18). This amounts to calculate the free energy of a domain wall separating 
the super phase and the normal phase with the bulk of super phase and that of the normal 
phase held in thermal equilibrium under an external magnetic field of the critical strength, 



The simplest ansatz of the solution to the Ginzburg-Landau equations that minimizes the 
domain wall free energy consists of only nonzero V, parallel to the external magnetic field, 
and nonzero components $o and ^s, which maintains the maximum symmetry spared by 
the boundary condition. Upon introducing dimensionless quantities via 

s^5x, $0 + = (Pou, X = ^(*o - ^*8) = V2<Pov, V = ^Acos e. 

the corresponding Ginzburg-Landau equations become 

- A" + lA(2u'' + v'') = 0, 

- ^u" + {A^ - l)u + l{2u^ + v^)u + Ip{u^ - v^)u - 0, 

- ^v" + - 4)^; + l(u'' + 5v^)v - |p(m2 _ v^)v = 0, (20) 

subject to the boundary conditions 

M I— > 0, w I— > 0, A' I— > 1 at s I— > — oo 

14 1— > 1, V I— > 1, ^4' 1-^ at s, I— > oo (21) 
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where k = | is the Ginzburg-Landau parameter, p = ll+^l with p — — | for mean field 
approximation, and the prime denotes the derivative with respect to s. In case of an ordi- 
nary superconductors with one component of the order parameter, Ginzburg-Landau found 
analytically [26] that the domain wall energy vanishes at 

K = = -L ^ 0.707, (22) 

which was later clarified by Arikosov [27] as the demarcation between the type I supercon- 
ductivity (k < Kc) and the type II one {k. > Kc)- The equations (20) are more complicated 
and depend on two dimensionless parameters, k and p. Nevertheless a set of inequalities have 
been established following the variational arguments [25]. We have shown that the domain 
wall energy is a decreasing function of k. and an increasing function of p. Consequently, 

The numerical solution of the equations (21) shows that 

Kc = 0.589 (23) 

for p = — |. While the color-flavor locked component $o dominates the bulk of the conden- 
sate, the unlocked component, shows up near the domain wall. 

Various fluctuation effects of the Ginzburg-Landau theory have also been considered in 
the literature, see e.g. [25] 



IV. CONCLUDING REMARKS 

In this lecture, I have reviewed some novel properties of the color superconductivity. A 
systematical approach to determine various physical quantities at weak coupling has been 
sketched. While the accuracy of the perturbative formula for the transition temperature 
(13) and that for the Ginzburg-Landau coefficients (16) is uncertain when extrapolated to a 
realistic quark matter, it is instructive to substitute in the obscrvationally attainable quark 
density for an order of magnitude estimation of the quantities of interests. 
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Taking the generic chemical potential of in the core of a neutron star, — 400MeV, as 

2 

a bench mark, we have = |^ ~ 1 and the mixing angle 9 = —5.6° for Nc = Nf = 3 
and A = Aqcd — 200MeV. Using eq. (16) and (19), we find the thermodynamical critical 



the color superconductivity is type I if ksTc < 14MeV and type II if ksT^ > 14MeV. It 
follows from (13) that /c^Tc ~ 3.5MeV with one-gluon exchange approximation and the color 
superconductivity is of type I. Beyond the one-gluon exchange approximation, ksTc could 
be high enough to cross over to the type II region, but it is unlikely to be a strong type II 
system. The critical magnetic field, which is by orders of magnitude higher than the typical 
magnetic field in a neutron star, is too strong to lead to observational distinctions between 
the two types of the color-superconductivity. 



I am grateful to W. Brown, I. Giannakis and J. T. Liu for collaborations. This work 
is supported in part by US Department of Energy contract number DE-FG02 -91ER40651- 
TASKB. 




Furthermore, the criterion (23) implies that 
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